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SUMMARY: In this paper, the integrated continuum radio spectrum of super¬ 
nova remnant (SNR) W44 was analyzed up to 70 GHz, testing the different emission 
models that can be responsible for its particular shape. The observations by the 
Planck space telescope made possible to analyze the high frequency part of radio 
emission from SNRs. Although the quality of radio continuum spectrum (a high 
scatter of data points at same frequencies) prevents us to make definite conclusions, 
the possibility of spinning dust emission detection towards this remnant is empha¬ 
sized. In addition, a concave-down feature, due to synchrotron losses, can not be 
definitely dismissed by the present knowledge of the integrated radio continuum 
spectrum of this SNR. 
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1. INTRODUCTION 

Supernova remnant W44 (G34.7-0.4, 3C392) 
represents one of the most interesting SNRs. It is a 
member of the mixed-morphology (thermal compos¬ 
ite) class of SNRs, characterized by a bright non- 
thermal shell-like radio morphology and centrally 
concentrated thermal X-ray emission (Rho & Petre 
1998; Vink 2012). W44 is a middle-aged SNR (possi¬ 
bly around 20 000 yrs old), that has an asymmetric, 
quasi-elliptical shell morphology, presumably due to 
expansion in an inhomogeneous interstellar medium 
(Cardillo et al. 2014). It is about half a degree in 
size (35' x 27') at a distance of around 3 kpc, and it 
is located in a complex region of the inner Galactic 
plane rich in both thermal and non-thermal sources 
(Castelletti et al. 2007). In fact, this remnant is well 
placed in the W48 molecular cloud complex, a rich 
star-forming region. 

SNR W44 constitutes one of the few demon¬ 
strated cases of an SNR-molecular cloud interact¬ 


ing system (Castelletti et al. 2007, and references 
therein). It is also detected in y-rays. Cardillo et al. 
(2014) concluded that y-ray emission is most proba¬ 
bly caused by the neutral pion decay (hadronic sce¬ 
nario). They also determined that the average gas 
density of the regions emitting 100 MeV - 10 GeV 
y-rays is relatively high (250 — 300 cm -3 ). Yoshiike 
et al. (2013) confirmed that the W44 SNR environ¬ 
ment consists of both molecular and atomic hydro¬ 
gen and they concluded that molecular clouds, which 
are likely associated with this remnant, surround the 
radio shell. Their analysis also supports a hadronic 
origin of the y-rays. 

The global integrated radio continuum spec¬ 
tral index of this SNR is a = 0.37 (Green 2014), 
where spectral index is defined as S v oc v ~ a , and S„ 
is the flux density. This is significantly less than 0.5, 
the value that is predicted by the test-particle diffu¬ 
sive shock acceleration theory (DSA, Urosevic 2014, 
and references therein). On the other hand, such 
values for spectral index are common for majority of 
mixed-morphology SNRs (Onic 2013, and references 
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therein). Uchida et al. (2012) reported another com¬ 
mon feature for mixed-morphology SNRs: radiative 
recombination continua of highly ionized atoms in 
X-ray spectrum. They concluded that the spectrum 
is well reproduced by a thermal plasma in a recom¬ 
bining phase. 

The spatially resolved radio spectral index 
study revealed that the detected localized absorp¬ 
tion has a negligible influence on the total integrated 
flux, and thus has no measurable impact on the inte¬ 
grated continuum spectrum (Castelletti et al. 2007). 
The most likely explanation of the particular local 
spectral inversion is the low-frequency radio contin¬ 
uum free-free absorption from ionized gas in the post 
shock region at the SNR/molecular cloud interface. 
An alternative explanation is that the thermal ab¬ 
sorption is occurring inside the boundaries of the co¬ 
incident H II region, along its periphery where the 
thermal electron density might be the highest. The 
spectral inversion is probably produced by a com¬ 
bination of both of these effects (Castelletti et al. 
2007). 

It is worth mentioning that there is no evi¬ 
dence in the radio continuum spectrum of any cou¬ 
pling between the associated pulsar PSR B1853+01 
and the surrounding SNR shell that could, for ex¬ 
ample, be observed as a gradual steepening from the 
pulsar to the shell (Castelletti et al. 2007). In ad¬ 
dition, the pulsar wind powers a small synchrotron 
nebula observed at radio frequencies and X-rays (An- 
derl et al. 2014, and references therein). 

A very interesting discovery of hard X-ray 
emission from SNR, not connected with associated 
pulsar or correspondent pulsar wind nebula, was re¬ 
ported by Uchida et al. (2012). They found that hard 
X-rays have an arc-like structure spatially-correlated 
with a radio continuum filament. Uchida et al. 
(2012) noted that the surface brightness distribution 
shows a clear anti-correlation with 12 CO( J = 2 — 1) 
emission from a molecular cloud. Finally, they con¬ 
cluded that the hard X-rays are most likely due to a 
synchrotron enhancement in the vicinity of the cloud. 
In fact, the localized non-thermal X-ray emission, as 
it is seen in the case of W44, probably reflects the 
physics of shock propagation in a clumpy medium. 
In an inter-clump gas with a lower density, cloud 
shocks can propagate much faster than in the dense 
clumps, so they can possibly accelerate high-energy 
electrons capable enough to emit X-rays (Lee et al. 
2015). 

Using Spitzer mid-IR observations, an H II re¬ 
gion, located just outside the southeastern limb of 
SNR W44, was detected by Castelletti et al. (2007). 
They showed that the combination of 8 and 24 fim 
Spitzer images reveals the hot dust grains in the 
Stromgren sphere, limited to the east by an annular 
photo dissociation region (PDR) dominated by poly¬ 
cyclic aromatic hydrocarbons (PAHs) emitting near 
8 /ini. Castelletti et al. (2007) noted that exactly at 
the interface between observed H II region and W44, 


as seen in the plane of the sky, a young stellar object 
is present. In fact, two massive young stellar objects 
are identified at the border of the Hii region that 
is evolving within a molecular cloud shocked by the 
SNR W44 (Paron et al. 2009; Ortega et al. 2010). 

Recently, the observations from a microwave 
survey of Galactic SNRs made by Plane /£] were 
published (Planck Collaboration Int. XXXI 2014). 
Planck observed the sky in nine frequency bands cov¬ 
ering 30 - 857 GHz with high sensitivity and a range 
of angular resolutions from 31' to 5' (Planck Collab¬ 
oration Int. XXXI 2014 and references therein). The 
Low Frequency Instrument covers the 30, 44, and 
70 GHz bands while the High Frequency Instrument 
covers the 100, 143, 217, 353, 545, and 857 GHz. 

Planck Collaboration Int. XXXI (2014) re¬ 
ported that the synchrotron emission from SNR W44 
is detected at levels above the emission from nearby 
unrelated regions at 30-70 GHz. They also noted 
that the measured 70 GHz flux density from Planck 
is somewhat lower than that expected from the radio 
power law, while the 30 GHz flux density is higher. 
Planck Collaboration Int. XXXI (2014), on the other 
hand, discounted the 30 GHz flux density due to pos¬ 
sible confusion with unrelated large-scale emission 
from the Galactic plane. Also, it must be noted that 
the Planck’s flux density measurements above 100 
GHz are contaminated by an unrelated foreground 
emission. Planck Collaboration Int. XXXI (2014) 
identified a particular structure at the eastern border 
of the radio SNR as a compact HII region, unrelated 
to the W44 SNR, but possibly a member of the same 
OB association as the progenitor. 

In this paper, we have analyzed the radio spec¬ 
trum of the W44 SNR up to 70 GHz, testing different 
emission models in order to replicate its particular 
spectral shape. 


2. ANALYSIS AND RESULTS 


Flux densities at different frequencies for SNR 
W44 were taken from Table 2 of Castelletti et al. 
(2007), Table 1 of Sun et al. (2011) and Table 3 of 
Planck Collaboration Int. XXXI (2014) in the fre¬ 
quency range from 610 MHz to 70 GHz. Only flux 
densities with errors < 20% were used in the anal¬ 
ysis and only data corrected to the scale of Baars 
et al. (1977) are taken from Castelletti et al. (2007). 
Finally, the misprint in Table 2 of Castelletti et al. 
(2007) regarding the data uncertainties taken from 
Altenhoff et al. (1970) is taken into account - uncer¬ 
tainties are set to 10% of the selected flux densities. 

The significant scatter of flux densities mea¬ 
sured at the same frequencies makes the integrated 
continuum radio spectrum difficult to analyze (see 
Figure 6 from Castelletti et al. 2007 and Figures in 
this paper). The data come from a wide variety of 
telescopes with different beam sizes, which can have 


1 Planck is a project of the European Space Agency (ESA) with instruments provided by two scientific consortia funded by 
ESA member states, with contributions from NASA (USA) and telescope reflectors provided by a collaboration between ESA 
and a scientific consortium led and funded by Denmark. 
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a significant impact on flux density measurements, 
since W44 lies within a particularly complex region 
close to the Galactic plane. The effects of possible 
contamination by the unrelated sources can be sig¬ 
nificant enough to account for the scatter in the radio 
continuum spectrum of W44. 

2.1. The synchrotron emission from the SNR 
W44 

The usual interpretation of the radio contin¬ 
uum spectrum of SNRs is that of a simple power law 
that arises from the synchrotron emission of charges 
accelerated by a diffusive shock acceleration (DSA) 
mechanism (Urosevic 2014, and references therein). 
In that case, the flux density is given by the following 
relation (Eq. 1) 

*S'[.Iy]( I/ ) = S[Jy](lGHz) P[ GHz ], (1) 

where a is the radio spectral index. On the other 
hand, the high frequency part of the W44 radio spec¬ 
trum clearly indicates that this simple model does 
not adequately explain the integrated continuum of 
this SNR. The lower value of the flux density at 70 
GHz could be a possible indication of high frequency 
spectral bending (Planck Collaboration Int. XXXI 
2014) that can be represented by (Eq. 2) 

-VlW = S[j y] (lGHz) p [ g “ z] e~"o, (2) 

where vq is a characteristic cut-off frequency. 

We performed least-squares fits to the radio 
spectrum (Fig. 1) using a simple synchrotron model 
(solid line) as well as a synchrotron model with expo¬ 
nential cutoff (dotted line). The thick and thin lines 
in Figure 1 correspond to fits without and with the 
inclusion of the Planck 30 GHz flux density, respec¬ 
tively. Planck data is represented with a diamond 
symbol, except for the data point at 30 GHz, which 
is labeled with a triangle. The best fitting param¬ 
eters for these two models are presented in Table 
1. The asterisk in Table 1 indicates the cases when 
Planck data at 30 GHz is included in the fit. 

It is known that the radio spectra of evolved 
SNRs could appear in a concave-down form. This 
kind of spectrum can be explained using DSA the¬ 
ory with the effect of synchrotron losses within the 
finite emission region. Generally, the concave-down 
form of the radio spectra should correspond to very 
old SNRs, for which electrons have had enough 
time to lose a significant amount of energy at the 
highest radio frequencies, and primarily to distant 
(mainly extragalactic) SNRs for which the limita¬ 
tion in telescope resolution leads to the capturing of 
radio emission from the sample of ” exhausted” elec¬ 
trons (Urosevic 2014, and references therein). The 
angular resolution of Planck ’s Low Frequency Instru¬ 
ment (LFI) at 30, 44 and 70 GHz is 33', 24', 14', 
respectively (Mandolesi et al. 2010). For a compar¬ 
ison, the angular resolution of Wilkinson Microwave 


Anisotropy Probe (WMAP) at 23, 33, 41, 61 and 94 
GHz is around 53', 40', 31', 21' and 13', respectively 
(Bennett et al. 2003). The resolution of WMAP at its 
lowest frequencies makes it unsuitable for our analy¬ 
sis of the W44.SNR. 

In addition, it must be noted that a few exam¬ 
ples of concave-down radio spectra of Galactic SNRs 
can be found in the recent literature, e.g. cases of 
SNRs S147 and HB21. Xiao et al. (2008) identified 
a spectral break at 1.5 GHz for SNR S147 and Pi- 
vato et al. (2013) reported steepening at around 6 
GHz in the case of SNR HB21. Finally, Planck Col¬ 
laboration Int. XXXI (2014) reported several other 
Galactic SNRs with possible high frequency spectral 
bending (such as IC443 and Puppis A). The high 
frequency spectral bending starts at relatively lower 
frequencies (several GHz) for S147 and HB21 than 
for SNRs W44, IC443 and Puppis A (several tens of 
GHz). 



Fig. 1. The weighted least-squares fit to the data for 
the synchrotron power law model (solid line) as well 
as for synchrotron power law m.odel with high fre¬ 
quency exponential cut-off (dotted line). Thick and 
thin lines correspond to the fits without and with in¬ 
clusion of Planck data at 30 GHz, respectively. Dia¬ 
mond symbols indicate Planck data; the point at 30 
GHz is shown as a triangle. 


2.2. The analysis of possible spinning dust 
emission towards the SNR W44 

The spinning dust emission is currently one of 
the most probable proposed mechanisms to explain 
the so called anomalous microwave emissior0 (AME) 
from the diffuse interstellar medium (ISM) of Milky 
Way (Planck Collaboration Int. XV 2014, and refer¬ 
ences therein). In fact, it is believed that the anoma¬ 
lous component of diffuse Galactic background is 
produced by electric dipole rotational emission from 
very small dust grains (Erickson 1957; Draine & 


2 Dust-correlated emission observed between around 10 — 100 GHz that can not be accounted for by extrapolating the thermal 
dust emission to low frequencies. 
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Lazarian 1998a). Small dust grains are likely to have 
a nonzero electric dipole moment due to the intrinsic 
dipole moment of molecules within the grain and un¬ 
even charge distribution. These grains will spin due 
to interaction with the ambient interstellar medium 
and radiation field, and thus radiate electromagnetic 
waves due to the rotation of their electric dipole mo¬ 
ment (Ali-Haimoud et al. 2009). To calculate the 
detailed frequency spectrum of spinning dust, one 
needs to integrate the emission over a distribution 
of grain sizes, electric dipole moments, and angular 
velocities (see Dickinson et al. 2014, and references 
therein for more details). A lot of recent papers deal 
with additional improvements of the basic emission 
model (Ali-Haimoud et al. 2009; Ysard et al. 2010; 
Ysard & Verstraete 2010; Hoang et al. 2010, 2011; 
Ysard et al. 2011). 

Spinning dust emission can generally con¬ 
tribute significantly at high radio continuum fre¬ 
quencies, especially around 10 - 100 GHz (Draine & 
Lazarian 1998ab; Ali-Haimoud et al. 2009; Steven¬ 
son 2014; Planck Collaboration Int. XV 2014). It 
is responsible for a characteristic bump in the high 
frequency part of the radio continuum. 

Scaife et al. (2007) fitted the radio spectrum of 
Galactic SNR 3C396 with a Warm Neutral Medium 
(WNM) spinning dust emission model (Draine & 
Lazarian 1998b), asserting for the first time the 
possibility of spinning dust (rapidly rotating dust 
grains) emission from the vicinity of an SNR. Al¬ 
though the claims of possible contribution of thermal 
bremsstrahlung at high frequencies were addressed 
by Onic et al. (2012) for SNR 3C396, the same is not 
applicable in the case of the particular shape of radio 
spectrum of W44. 

In their study of AME in Galactic clouds, 
Planck Collaboration Int. XV (2014) listed W48 
complex (that contains W44) as one of their can¬ 
didate regions. However, their conclusions, based on 
a significance level of the AME detection, suggest 
that the excess emission, for this object, is not sta¬ 
tistically highly significant. On the other hand, they 
emphasized that many of the so called statistically 
semi-significant AME regions (those at 2 — 5cr), like 
W48, are likely to be real detections of AME. Planck 
Collaboration Int. XV (2014) also noted that in gen¬ 
eral, the emerging picture is that the bulk of the 
AME (most probably due to spinning dust) is com¬ 
ing from the polycyclic aromatic hydrocarbons and 
small dust grains from the colder neutral interstel¬ 
lar medium phase. They also suggested that AME 
comes from the molecular cloud dust or PDR, but 
not from H II regions themselves. That is consistent 
with a general belief that PAHs are destroyed in H II 
regions, and the AME emissivity is lower in the ion¬ 
ized phase of the interstellar medium. 

However, Hensley & Draine (2015) voiced 
their concern on the validity of claims about 
AME/PAHs association. They proposed that one 
possibility is that AME could be in fact a spinning 
dust emission that arises primarily from very small 


grains that are not PAHs. Of course, AME could be a 
combination of the emission from spinning PAHs and 
non-PAH spinning dust, as well as the thermal dust 
emission, such as magnetic dipole emission (Hensley 
& Draine 2015, and references therein). Hensley & 
Draine (2015) did not find the evidence for strong 
AME correlation with free-free or CO emission, but 
they emphasize the plausibility of correlation with 
synchrotron emission. If the AME arises from spin¬ 
ning ultra small grains, it might be enhanced in SNRs 
in which grains are violently shattered. In spite 
of that, it can also be suppressed as the result of 
the destruction of very small grains by sputtering 
(Lakicevic et al. 2015). 

In a recent paper, Irfan et al. (2015) analyzed 
the region that encompasses SNR W44. The inner 
aperture size used for their aperture photometry was 
set to 60' and centered around the Galactic coordi¬ 
nates, l = 34? 8, b = — 0?5. Irfan et al. (2015) 
demonstrated the presence of AME associated with 
that region. They also ruled out the possibility that 
the excess emission near 30 GHz was from a nearby 
ultra-compact H n region. 

As the proper physical conditions may ex¬ 
ist, it is worth checking if spinning dust emission 
is significant enough to shape the radio spectrum of 
W44 SNR near 30 GHz. To that end, the SpDust 
code, ver. 2.01 (Ali-Hahnoud et al. 2009; Silsbee et 
al. 2011) was used. The radio continuum spectrum 
from 610 MHz to 70 GHz was fitted by the sum of 
synchrotron radiation represented by the power law 
with spectral index a and the spinning dust emission 
S sp d(v,n,T) (Eq. 3) 

S[jy] M = ^(lGHz) ^-“ z] + S spd (u- n, T). (3) 

A spinning dust emission spectrum from a 
WNM model (assume nn = 0.4 cm -3 , T = 6000 K) 
was used in fitting the overall radio continuum spec¬ 
trum, as well as several other standard types of envi¬ 
ronments that were used for comparison: Warm Ion¬ 
ized Medium (WIM, «h = 0.1 cm -3 , T = 8000 K), 
Cold Neutral Medium (CNM, ?rn = 30 cm -3 , T = 
100 K), and Molecular Cloud (MC, «h = 300 cm -3 , 
T = 20 K). These idealized ISM phases are defined in 
Table 1 of Draine & Lazarian (1998b). The param¬ 
eters for the grain size distribution are taken from 
Table 1 of Weingartner & Draine (2001) in accor¬ 
dance with the work of Ali-Haimoud et al. (2009) - 
see Figure 14 of their work. Due to the low qual¬ 
ity of radio spectrum we must confine ourselves to 
the qualitative analysis using standard ISM environ¬ 
ments. 

The weighted least-squares fit is calculated us¬ 
ing the MPFIT0 (Markwardt 2009) package written 
in IDL for all of the fits presented in this paper, with 
starting values estimated from the data. MPFIT also 
provides estimates of the ler uncertainties for each 
parameter, taken as the square root of the diagonal 
elements of the parameter covariance matrix. We 
must note here that while modeling spinning dust 
emission can yield very good fits, one must exercise 


3 http://purl.com/net/mpfit 
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caution in interpreting these results since the fitted 
parameters can become unphysical. 

Also, it is worth mentioning that in the case 
of non-linear model fitting, the number of degrees 
of freedom is generally unknown, i.e. it is not possi¬ 
ble to compute the value of reduced % 2 , or adjusted 



« [GHz] 



« [GHz] 


R 2 (Andrae et al. 2010). In non-linear models, the 
parameter k = N — p (where N is the number of 
data points and p is the number of model parame¬ 
ters) does not always represent the exact number of 
degrees of freedom (Andrae et al. 2010). 



■z [GHz] 



« [GHz] 


Fig. 2. The weighted least-squares fit to the data for the synchrotron power law model with inclusion of 
spinning dust emission (Eq. 3) is presented on the upper graphs and the weighted least-squares fit to the 
data for the synchrotron power law model with high frequency exponential cut-off with inclusion of spinning 
dust emission (Eq. f) is presented on the lower graphs. The left and right graphs correspond to the data 
samples with flux density at 30 GHz excluded and included, respectively. Diamond symbols indicate Planck 
data and triangle symbol represents the data point at 30 GHz. Different lines correspond to different models 
of idealized phases of interstellar medium: WIM - solid line, WNM - dotted line, CNM - dashed line, MC - 
dashed dot line. 


The results of weighted least-squares fit to 
data by Eq. 3 are presented in left and right up¬ 
per graphs of Figure 2. The left and right graphs 
correspond to the data samples with flux density at 
30 GHz excluded and included, respectively. Dia¬ 
mond symbols indicate Planck data and the triangle 
symbol represents data at 30 GHz. Different lines 
correspond to different models of idealized phases of 
interstellar medium: WIM - solid line, WNM - dot¬ 
ted line, CNM - dashed line, MC - dashed dot line. 


The best fitting parameters as well as corresponding 
X 2 ( k ) for different models are presented in Table 1 
(star symbol indicates inclusion of the 30 GHz Planck 
data in the analysis). 

The large scatter of flux densities measured at 
the same frequencies prevents us from making a firm 
discrimination between the models we are consider¬ 
ing. We cannot rule out the significance of spinning 
dust emission from this SNR. A model that involves 
the WNM condition is statistically the most probable 
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(supported by the standard analysis of y 2 /fc values) 
and also physically plausible for the particular SNR 
environment. In addition, this model (dotted line in 
Figure 2) passes through the error bars of the data 
point at 30 GHz even when that point is excluded 
from the analysis. As it was noted by Planck Col¬ 
laboration Int. XXXI (2014), this point is possibly 
contaminated by unrelated large-scale emission from 
the Galactic plane. Although such a claim can not 
be fully dismissed, the strong possibility exists that 


the spinning dust emission is in fact, at least par¬ 
tially, responsible for the apparent bump around 30 
GHz. 

For completeness, the spinning dust flux den¬ 
sity estimate 5 sp d(30), as well as its ratio to the to¬ 
tal flux density /(30) = S' sp d(30)/S'(30) at 30 GHz 
is given in Table 1 for all considered models. The 
spinning dust fractions /( 30) are in accordance with 
the results obtained for other AME sources (~ 0.5) 
in Planck Collaboration Int. XV (2014). 


Table 1. The best fitting parameters for different models. 


Equation 

a 

i/ 0 [GHz] 

5 spd (30) [Jy] 

/(30) 

X^ (k) 

1 

0.433 ±0.014 

- 

- 

- 

268.18 (26) 

1* 

0.409 ±0.013 

- 

- 

- 

322.11 (27) 

2 

0.388 ±0.024 

234 ± 109 

— 

— 

262.84 (25) 

2* 

0.358 ±0.023 

220 ±91 

- 

- 

315.13 (26) 

3 (WIM) 

0.467 ±0.018 

— 

51.6 ± 12.9 

0.50 ±0.20 

251.76 (25) 

3* (WIM) 

0.473 ±0.017 

— 

63.5 ± 7.7 

0.55 ±0.12 

253.13 (26) 

3 (WNM) 

0.473 ±0.018 

— 

78.7 ± 18.0 

0.60 ±0.24 

248.63 (25) 

3* (WNM) 

0.470 ±0.017 

- 

70.6 ±8.3 

0.57 ±0.12 

248.90 (26) 

3 (CNM) 

0.465 ±0.018 

— 

40.1 ±10.6 

0.43 ±0.18 

253.60 (25) 

3*(CNM) 

0.473 ±0.017 

— 

57.1 ±7.2 

0.52 ±0.12 

258.53 (26) 

3 (MC) 

0.456 ±0.019 

— 

9.7 ±5.0 

0.15 ±0.10 

264.09 (25) 

3*(MC) 

0.460 ±0.019 

- 

20.3 ±4.7 

0.27 ±0.10 

303.06 (26) 

4 (WIM) 

0.408 ±0.025 

136 ± 46 

63.2 ± 13.7 

0.55 ±0.23 

241.14 (24) 

4* (WIM) 

0.409 ±0.025 

130 ± 42 

68.4 ±8.0 

0.58 ±0.15 

241.35 (25) 

4 (WNM) 

0.428 ±0.027 

196 ± 88 

81.7 ± 18.3 

0.62 ±0.26 

242.93 (24) 

4* (WNM) 

0.424 ±0.026 

203 ± 91 

71.3 ±8.4 

0.58 ±0.15 

243.34 (25) 

4 (CNM) 

0.403 ±0.025 

127 ±42 

52.0 ± 11.5 

0.51 ±0.21 

241.83 (24) 

4* (CNM) 

0.402 ± 0.025 

112 ±33 

64.0 ± 7.7 

0.56 ±0.15 

243.75 (25) 

4 (MC) 

0.374 ±0.028 

61 ±22 

31.5 ±7.9 

0.42 ± 0.22 

243.32 (24) 

4*(MC) 

0.354 ±0.031 

47 ± 17 

43.4 ±8.1 

0.52 ±0.23 

268.05 (25) 


*For the case when Planck's data point at 30 GHz is included in the analysis. 


Additional support for a spinning dust emis¬ 
sion hypothesis can be found in its correlation with 
thermal dust emission, especially at IRAS (Infrared 
Astronomical Satellite) wavelengths. Planck Collab¬ 
oration Int. XV (2014) found that AME sources gen¬ 
erally have a 12 /mr/25 /im ratio ss (0.6 — 1.0). This 
is interpreted as a confirmation of the spinning dust 
model, where the very small grains are responsible 
for the bulk of the AME. Although the W44 SNR was 
not detected by Saken et al. (1992), Arendt (1989) 
found a 12 /im/25 /mi ratio of ss 1. However, pos¬ 
sibility of ionic line contamination of the particular 
IRAS bands can not be fully dismissed (Oliva et al. 
1999; Reach & Rho 1996). 

The comparison was also made between the 
spinning dust flux density estimates at 30 GHz 
5s P d(30) for different models to the 100 /mi (3000 
GHz) flux density 5(3000) as given in Arendt (1989). 


The ratio 5 sp d(30)/5(3000) spans a large range of 
values for the different environmental models: (7 — 
60) x 10~ 4 . Bearing in mind the roughness of our 
analysis, as well as high uncertainty of the IRAS 
100 /im flux density (around 50%), our values for 
5 sp d(30)/5(3000) are in a rough agreement (slightly 
higher) with other determinations for AME sources 
((1 — 15) x 10 -4 , Planck Collaboration Int. XV 2014; 
Hensley & Draine 2015, and references therein). 

Finally, it is tempting to check whether the fits 
will improve if the model that incorporates the sum 
of power law with exponential cut-off and spinning 
dust emission is used (Eq. 4) 

S[jy]M =S^ c (lGHz) u^ z] e-^-o + 5 spd (^;n,T). 

( 4 ) 

The lower graphs of Figure 2 represent the 
weighted least-squares fit to the data for the syn- 
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chrotron power law model with high frequency expo¬ 
nential cut-off with inclusion of spinning dust emis¬ 
sion (Eq. 4). The left and right graphs correspond 
to the data samples with flux density at 30 GHz ex¬ 
cluded and included, respectively. The best fitting 
parameters as well as corresponding x 1 2 (fc) for dif¬ 
ferent models are given in Table 1. 

The average x 2 /fc is slightly less for Eq. 4 than 
for Eq. 3 (see Table 1). On the other hand, statis¬ 
tical discrimination between different environmental 
conditions used for Eq. 4 is in this case even more 
inconclusive. Still, detection of hard X-ray emis¬ 
sion from W44, not connected with associated pul¬ 
sar or its nebula, reported by Uchida et al. (2012), 
can cast doubt that high frequency synchrotron cut¬ 
off at around 100 - 200 GHz (see Table 1) is gen¬ 
uine and is primarily responsible for a lower flux 
density at 70 GHz. On the other hand, Lee et al. 
(2015) recently explored the non-thermal emission 
mechanisms of dynamically evolved SNRs. They ex¬ 
plored two scenarios of particle acceleration, either 
a re-acceleration of Galactic cosmic rays, or an effi¬ 
cient nonlinear diffusive shock acceleration of ther¬ 
mally injected particles by a fast radiative cloud 
shock. Lee et al. (2015) emphasized that if suffi¬ 
ciently strong magnetic turbulence is present in the 
molecular cloud, the re-acceleration scenario of the 
non-thermal emission agrees well with the broad¬ 
band spectrum of SNR W44 and that the apparent 
discrepancy at 30 GHz can possibly arise from the 
anomalous microwave emission from small spinning 
dust grains. Their predicted radio spectrum shows a 
spectral softening above around 10 GHz due to syn¬ 
chrotron loss, and is consistent with the 70 GHz flux 
density from Planck. 

One must bear in mind that SNR W44 is actu¬ 
ally placed in a very complex region (such that phys¬ 
ical conditions inside the remnant may significantly 
vary) so that the more detailed models than these 
presented in this paper should be more adequate. 
In fact, thorough radiative transfer modeling is re¬ 
quired to further our understanding of the spinning 
dust emission and its significance. Finally, statisti¬ 
cally insufficient number of data prevents us to make 
a firm quantitative discrimination between different 
models. In addition to the general need for more 
data, the much more improved resolution (< 1') of 
new observations at radio continuum frequencies be¬ 
tween 10 and 100 GHz would also be beneficial. 


3. CONCLUSIONS 


In this paper, the integrated radio spectrum 
of SNR W44 was analyzed up to 70 GHz testing the 
different emission models that can be responsible for 
its particular shape. The main conclusions are: 

(1) Planck observations made it possible to ana¬ 
lyze the high frequency part of radio emission 
from SNRs. 

(2) Although the quality of radio continuum spec¬ 
trum (a high scatter of data points at same fre¬ 
quencies) prevents us to make definite conclu¬ 


sions, the possibility of spinning dust emission 
detection towards this remnant is emphasized. 
In fact, the spinning dust emission is proposed 
to be, at least partially, responsible for the ap¬ 
parent bump in the radio continuum around 
30 GHz. 

(3) In addition, a concave-down feature in the 
radio spectrum, due to synchrotron losses, 
can not be definitely dismissed by the present 
knowledge of the radio continuum spectrum of 
this SNR. 

(4) Finally, the more data at different frequencies 
between 10 and 100 GHz at much better res¬ 
olution than Planck are needed to make firm 
conclusions about the contribution of partic¬ 
ular radiation mechanisms responsible for an 
observed shape of the radio spectrum of SNR 
W44. 
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ON THE RADIO SPECTRUM OF SNR W44 


O MHTErPAJIHOM HEIIPEKMUHOM PAHMO-CIIEKTPy OCTATKA 
CynEPHOBE W44 (G34.7-0.4): HOBA CA3HAEBA nOMOTiy I1JIAHKA 
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y obom pa^y je aHajm3HpaH HHTerpajimi 
HenpeKna;HH paario-cneKTap ocTaTKa cynepHOBe 
(OCH) W44 CBe p,o 70 GHz m TecTupami 

Cy pa3JIH^HTH eMHCHOHH MO^eJIH Kojn Mory 

y3poKOBaTH meroB CBojeBpcHH o6jihk. noc- 
MaTpaita peanH30BaHa cbcmhpckhm TejiecKonoM 
njiaHK OMorytiHJia cy .zja ce aHajnoupa bh- 
coKO())peKBeHTHn &eo panno-eMHcnje OCH. Haico 
je KBajiHTeT HenpeKH^Hor pa^HO-cneKipa (Be- 


jihko pacHnaae BpeAHOCTH Ha hcthm (j) p eKB eHH,H- 
jaMa) TaKaB n& HaM OHeMoryHaBa H3Be,neMO 
HBpcTe 3aKJi>yHKe, HCTaKHyTa je HajBepoBaT- 
Hnja a;eTeKpnja eMHCHje poTHpajyfie npammie y 
npaBpy OBor ocTaTKa. y3 to, TpeHyTHO no3Ha- 
Baae HHTerpajiHor HenpeKH,nHor panno-cneKTpa 
obot OCH Hnje ^obojbho na 6m ce .hc^hhhthbho 
M orjia ocnopHTM TBp^aa o KOHKaBHOM KpHBJtemy 
cneKTpa ycjien chhxpotpohckhx ry6HTafca. 
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